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Q Abstract —Distinction of OFDM signals from single carrier 
signals is highly important for adaptive receiver algorithms and 
signal identification applications. OFDM signals exhibit Gaussian 
characteristics in time domain and fourth order cumulants of 
Gaussian distributed signals vanish in contrary to the cumulants 
of other signals. Thus fourth order cumulants can be utilized for 
OFDM signal identification. In this paper, first, formulations of 
the estimates of the fourth order cumulants for OFDM signals 
are provided. Then it is shown these estimates are affected 
significantly from the wireless channel impairments, frequency 
offset, phase offset and sampling mismatch. To overcome these 
problems, a general chi-square constant false alarm rate Gaus- 
sianity test which employs estimates of cumulants and their 
covariances is adapted to the specific case of wireless OFDM 
signals. Estimation of the covariance matrix of the fourth order 
cumulants are greatly simplified peculiar to the OFDM signals. 
A measurement setup is developed to analyze the performance 
of the identification method and for comparison purposes. A 
parametric measurement analysis is provided depending on 
modulation order, signal to noise ratio, number of symbols, and 
degree of freedom of the underlying test. The proposed method 
outperforms statistical tests which are based on fixed thresholds 
or empirical values, while a priori information requirement and 
complexity of the proposed method are lower than the coherent 
identification techniques. 


I. Introduction 

Modulation identification has always been a part of commu¬ 
nications surveillance systems and electronic warfare, however 
gained more importance with the adaptive communication 
systems. Seamless and efficient communications should be 
achieved between a wide variety of communications systems 
which may employ different modulation techniques. Therefore, 
simple, robust, and fast modulation identification methods 
are required at the implementation of adaptive and intelli¬ 
gent receivers. Blind receiver design algorithms and signal 
transmission procedures can also benefit from the modulation 
recognition information. 

Improving the efficiency of the wireless communications 
systems, on the other hand, has substantial importance due to 
the increasing demand for higher capacity. Ascending number 
of users, networks, and deployment of a wide variety of 
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systems led the wireless communications towards heteroge¬ 
neous networks, small cell deployments and implementation 
of dynamic spectrum access techniques. Therefore emerging 
technologies require as much information as possible about 
the communications medium. Wireless communications mea¬ 
surement systems should also adapt themselves to the new 
wireless communications paradigm by providing more infor¬ 
mation such as modulation type, carrier number, bandwidth, 
and other waveform parameters automatically. To this end, 
radiocommunication sector of International Telecommunica¬ 
tion Union (ITU-R) initialized an incentive, and published the 
“Technical identification of digital signals” document as the 
main guideline for the wireless communications measurement 
systems to identify unexpected or unlicensed emissions, to de¬ 
tect interference sources, and to provide waveform parameters 
via signal identification H]. 

A. Related Work 

When the contemporary communications systems are con¬ 
sidered, orthogonal frequency-division multiplexing (OFDM) 
is a multi-carrier (MC) multiplexing scheme which maintains 
adaptive communications features by employing sub-carriers 
in a flexible way. Deployments of OFDM based systems are 
increasing rapidly. Distinction of OFDM signals from the 
single carrier (SC) signals is an important issue from the 
adaptive transceiver design, signal identification and spectrum 
monitoring perspectives. Initial research on this field indi¬ 
cates that modulation identification or recognition methods 
for OFDM systems can be classified into two groups; co¬ 
herent techniques which include maximum likelihood (ML) 
modulation classifiers along with other methods and non¬ 
coherent statistical methods which do not require a priori 
information, in general. Coherent techniques are based on the 
assumption that certain characteristics of the OFDM signals 
are known to the receiver. Statistical methods mostly depend 
on high order statistics (HOS) of cumulant or moment esti¬ 
mates of the received signals. The background of cumulant 
based modulation classification laid down by ||2l for fourth 
and sixth order cumulants. Later, fourth order cumulants are 
adapted to OFDM modulation identification Ei-iia . Moreover, 
fourth order moments of OFDM signals are combined with 
the amplitude estimation to improve signal detection in 0. 
On the other hand, support vector machines in combination 
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with specific fourth order cumulants are used for the same 
purpose in Q. Tabular values of high order moments which 
were implemented for SC modulation identification in ||8l are 
adapted to a minimum mean-square error (MMSE) estimator 
for OFDM signals in M- A low complexity single carrier 
modulation identification algorithm based on fourth order 
cumulants for multiple-input and multiple-output systems is 
also proposed in IfTOl . 

The ML modulation classification algorithms are known 
with better estimation results when compared to statistical 
methods, however they suffer from computational complexity 
and require channel estimation IfTTI . On the other hand, boot¬ 
strap technique is introduced as a pattern recognition based 
method in ifT^ . This method conducts SC modulation separa¬ 
tion along with the SC/MC distinction. A wavelet transform 
algorithm is proposed to extract transient characteristics of 
OFDM signals in IfTSlI . Moreover, time-frequency representa¬ 
tion of the wireless signals is employed to distinguish OFDM 
signals from SC frequency-shift keying (FSK) and phase-shift 
keying (PSK) signals in lfT4l . OFDM signals have unique 
second order cyclostationary features and a feature extraction 
algorithm based on the magnitude of the cyclic cumulants 
of OFDM signals is implemented in US). Furthermore, sec¬ 
ond order cyclostationary features of block transmitted-single 
carrier linearly digitally (BT-SCLD) modulated signals are 
utilized for distinction from OFDM signals and for parameter 
estimation in M- Another second order cyclostationary algo¬ 
rithm to detect and classify both worldwide interoperability for 
microwave access (WiMAX) and long term evolution (LTE) 
signals is developed in the context of cognitive radio systems 
and spectrum awareness concept in ini. An analysis and 
comparison of second order cyclostationary, matched filter, 
MMSE, and normalized kurtosis methods can be found in ifTSl 
and for a general overview of coherent techniques, m can 
be consulted. 

OFDM signals also exhibit time-domain Gaussianity. This is 
due to fact that the assignment of random data over orthogonal 
sub-carriers in a simultaneous way can be assumed as a com¬ 
position of large number of independent, identically distributed 
(i.i.d.) random variables and central limit theorem implies 
Gaussian distribution for large enough data sets. A Gaussianity 
test based on empirical distribution function is introduced for 
OFDM systems in IMI . In addition to that, a general chi- 
squared constant false alarm rate (CFAR) Gaussianity test 
based on the estimates of the third and fourth order cumulants 
m is adapted to OFDM signals without considering the 
channel effects in Il22ll . 

B. Proposed Method 

When the coherent methods are taken into account, there 
are two main drawbacks when compared to the statistical 
methods: first, the dependence on the specific signal parame¬ 
ters such as symbol duration, cyclic-prefix size, active carrier 
numbers for the identification of each signal type and standard. 
Second, the algorithm complexity can increase the software 
and/or hardware complexity thus complicate the feasibility of 
implementation of the method. Moreover, complexity issues 


can lead to extension of identification duration. On the other 
hand, it is known that statistical methods are simpler than 
the coherent methods however they are highly susceptible to 
the changes in the wireless channel conditions along with 
other impairments such as frequency offset, phase offset 
due to sampling mismatch and other factors ll2^ . Moreover, 
the statistical tests in the literature depend on either fixed 
thresholds or empirical values of various HOS. However HOS 
can significantly fluctuate under wireless channel conditions 
and the identification performances of these tests degrade 
significantly. 

It is known that considering all combinations of possible 
lags is a computationally consuming process for cumulant 
based statistical classification methods. However peculiar to 
the OFDM signals, inclusion of certain lags in the analysis is 
sufficient to retain Gaussianity. Therefore in this study, first 
the fourth order cumulant equations for the OFDM signals 
are provided for limited number and combination of lags and 
lag distances. Second, in ll^ a general analysis of the effects 
of wireless impairments over the fourth order cumulants of 
OFDM signals were conducted. Furthermore in this paper, the 
fourth order cumulant estimates including the wireless effects 
are provided. These equations imply fundamental impact of 
impairments over the cumulant values. Third, an OFDM signal 
identification method which is based on a modified chi-squared 
Gaussianity test is proposed to distinguish MC OFDM signals 
from SC signals. The method employs the estimates of the 
fourth order cumulants and their covariance matrix in the 
decision process. Therefore, a certain level of protection is 
achieved against the effects of multipath fading channels, when 
compared to the statistical methods which are based on only 
the estimates of cumulants or moments of wireless signals. 
Moreover as detailed in Section |III1 instead of fixed numerical 
values, the proposed identification method utilizes the degree 
of freedom of the underlying chi-squared test which provides 
relative isolation from the effects of wireless impairments. 

Fourth, the proposed method includes the computation of 
the estimate of the covariance matrix of fourth order cumu¬ 
lants. In general, this is a complex and detailed operation 
because the combinations of products of the second and 
fourth order moments are included in the estimation process. 
However, in this paper the estimation process is reduced down 
to three particular equations peculiar to the OFDM signals. 
The only parameter required is the symbol duration which 
determine the degree of freedom of the modified chi-squared 
test. The numerical values that satisfy the confidence levels 
for correct detection can be found iteratively thus the test can 
become fully blind. Fifth, a measurement setup is developed 
and the proposed method is implemented to achieve parametric 
performance analysis in a realistic way. A signal generation 
software is developed and integrated to the signal generator 
to enable parametric flexibility during the signal transmission 
process. Such an approach led to detailed parametric analysis 
of the proposed method based on modulation order, signal 
to noise ratio (SNR), number of symbols, and degree of 
freedom of the underlying test against FSK, PSK, and quadra¬ 
ture amplitude modulation (QAM) SC signals. The detection 
performance of the proposed method is also compared with 
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three statistical methods which are based on HOS and a 
coherent identification method which employs second order 
cyclostationary features of OFDM signals for identification. 

The OFDM signal model which is used at the imple¬ 
mentation of the proposed method is defined in Section |II] 
The fourth order cumulant expressions for OFDM signals, 
discussion on the effects of wireless impairments on the 
cumulants, the modified chi-squared test which is employed 
by the proposed method, the decision mechanism, and the flow 
of the identification method are provided in Section |III1 The 
estimates of the covariances of the cumulants are derived in 
the Appendix. The measurement setup is detailed in Section 
lYl Parametric performance analysis of the proposed signal 
identification method based on the measurements is given in 
Section |V] Finally, conclusion constituted the Section Ivn 

II. Signal Model 

Inverse discrete Fourier transform (IDFT) is employed by 
OFDM systems for transmission. The baseband continuous¬ 
time OFDM signal is given by 

K 

Z ; 2na{k)i 

Sn(k)e^ -TG<t <Td, ( 1 ) 

where j at the exponential term is the imaginary unit, a is 
the set of K subcarriers in total, S„{k) is the n''^ OFDM data 
symbol which is transmitted at the k''' subcarrier from the 
set of a{k), To is the useful data duration, Tc is the cyclic 
prefix (CP) duration, and the total OFDM symbol duration is 
Ts - To + Tq. The transmitted total signal can be written as 

oa 

x{t) = ^ x„(t-nTs). (2) 

n=-oo 

The signal is modulated after digital analog (D/A) conver¬ 
sion and passed through mobile radio channel. The channel 
can be modeled as a time-variant linear filter 

L 

h(t) = hm{t)5{t - Tm), (3) 

m=l 

where L is the number of taps and t, is the excess delay 
for each tap. It is assumed that the taps are sample spaced 
and the channel is constant for a symbol but time-varying 
across multiple OFDM symbols. Therefore, baseband model 
of the received signal after down conversion also considering 
the channel model given in eqn. (O can be given by 

L 

y{t) = Y “ T,„)h,„(t) -i- n(f), (4) 

m=l 

where ^ is the carrier frequency offset due to inaccurate 
frequency synchronization, 0, is the phase offset, and n{t) 
corresponds to additive white Gaussian noise (AWGN) sample 
with zero mean and variance of cr^. The received signal is 
sampled with sampling time of Af at the analog to digital 
converter (ADC) and discrete-time received signal can be 
represented by 

L 

y{i) = _ T-^)/j^(/Af) -H n{i), (5) 

m=l 


where i - 1,..., A, Afl, and Af are normalized phase and 
carrier frequency offsets. Models of FSK, PSK, and QAM 
modulated SC signals are not listed here, but can be found 
in 13 and references therein. Finally, at the initial stage of 
the identification process, it is assumed that the receiver is 
focused on a single channel which may be utilized by a 
SC or a MC OFDM signal comprised of sub-carriers, the 
bandwidth and carrier frequency are estimated, the signal is 
filtered, down-converted and sampled beforehand. It is also 
assumed that in case of signals that are overlapping on wireless 
spectrum, before the initialization of the identification process, 
the proposed method can benefit from the signal separation 
processes introduced in 1241 . 

III. Proposed Identieication Method 

The proposed method benefits from the high order cumu¬ 
lants of the received signals for signal identification. The 
cumulants of a process are defined as the generalization of 
the autocorrelation function E{y(i)y(i -i- ii)}. Thus, the general 
formulation of third and fourth order cumulants for the sta¬ 
tionary random processes with zero mean is given by 

Csyih, h) - E{y(0y(i + h)y{i + h)] (6) 

and 

C4y{i\,h, h) - E{y(0y(i + i\)yii + h)y{i + I's)) 

- C2y{i\)C2y{i2 - h) - C2y(i2)C2y(h - ii) 

- C2y{h)C2y{i2 - h), (7) 

where C 2 y{i\) - E{y{i)y{i + I'l)). The third and fourth order cu¬ 
mulants have different characteristics; the A:th order cumulants 
will vanish for A: > 3 if y{i) is Gaussian l25l . However, third 
order cumulants can converge to zero although y{i) is non- 
Gaussian but symmetrically distributed l2^ . Therefore, third 
order cumulants will be ignored and fourth order cumulants 
will be analyzed for OEDM signal identification. On the other 
hand, when the calculation of the estimates of the cumulants 
are considered, the estimated cumulants, Cky, should be abso¬ 
lutely summable Il27l and the signal model for wireless OFDM 
signals which is given in eqn. © satisfies this condition. 
Under these conditions, the estimate autocorrelation becomes 
C 2 y{ii) - jj y(0y(i + h), and 4th-order cumulant estimate 
is given by 

j N-h 

C4y(iui2, h) - y(i)y(i + ii)y(i + u)y(i + h) 

1=1 

- C2y(il)C2y(i2 - h) - C2y(i2)C2y(h - il) 

- C2y{h)c2y{i2 - i\), {i\,i2, h) ^ 7^. (8) 

The region of all combinations of all lags is given by 
- {0 < 4-1 < ■ ■ ■ < i'l < oo). However, it is shown in Il28l 
that cumulant lags that comprise distribution characterization 
of the series are finite and it’s not required to estimate the 
cumulant values for the majority of the possible lag combina¬ 
tions. Therefore, the region of the lags that should be taken 
into consideration for the analysis is given by 

= {0 < is < h < i\ < N], 


(9) 
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and the sample estimates in eqn. ([8]l should only be calculated 
over as defined in eqn. (|9]l- The Gaussianity test will use 
lags of C 4 y that can be collected into an x 1 vector. This 
vector constitutes a 3 dimensional triangular region and its 
length is given by Nc - N(N + 1)(N + 2)/6. 


A. Fourth Order Cumulants for OFDM Signals 

When the wireless signals are considered, the length of the 
data can be orders of thousands or tens of thousands of sam¬ 
ples depending on the sampling rate and the recording time. 
Defining the lag region and computing the components of 
vector C 4 y can become impractical in many cases. 

The data set can be divided into smaller sections and 
processed in parallel or the data set can be shortened. However, 
even though the choice of N and consequently Nc is appli¬ 
cation dependent, in practice when a strong non-Gaussianity 
(e.g., due to the multipath fading channel) present in short 
records i.e, N < 200, this will weaken the Gaussianity 
assumption and will lead to unreliability for the signal identi¬ 
fication methods. Moreover, even for the shortest data sets, 
the size of the C 4 y vector can become very large and the 
computation of the cumulants can be impractical. However, 
Il29l and Il22l showed that, first, shrinking the lag region 
into = {/3 - 0 < i 2 — ii < M] does not lead to 
loss of significant distribution information for communications 
signals. The terms that are left over have very small influence 
on the features of the modulation type of the signal when 
compared to the set in the hand. Second, the lag distance 
can also be limited with M « l.STs where Ts is the symbol 
duration for the signal under test. Therefore the cumulants 
should be evaluated only for the set of lags of (A, A, 0) where 
ii - 12 - A e [0, l.Sr^]. Under these conditions eqn. ® can 
be written as 

N-A 

C 4 y(A, /1,0) = — ^ y{i)yii + A)y(i + A)y(i) 

^ 1=1 

- C 2 y( 0 )c 2 y( 0 ), (ii,i2,i3)e/f- (10) 

and the open form becomes 


C4y(A,A,0) 


. N~A r , N-A 

^ X - nTj 

?~i /=i 

1 ir 1 ^ 

-Y^y(i)y(i + A) - J^y(i)y(i-A) 


i = A+l 


l=l 


( 11 . 1 ) 

( 11 . 2 ) 

(11.3) 

( 11 ) 


When the signal model in the eqn. (|5]l is replaced into the 
eqn. (II 1.1b . the statements that are given in the eqn. (fThl l is 
achieved. Instead of expanding whole eqn. (|5]l, eqn. (fThl l is 
sufficient for the analysis at the next step. The results can 
then be generalized for the rest of the eqn. dD. Note that 
the noise component n{i) is independent from the rest of the 


terms that constitute the received signal and vanishes due 
to its Gaussianity. Therefore any term that included noise 
component is eliminated from the expressions given in eqn. 
(fThl) . When further mathematical operations on eqn. (fThl) 
is sought, the HOS theory ll2^ indicates that if rju, where 
u - 1,.. .,k- 1, are constants then 


Cky(T]iiu---,m-ik-i) = YXrjAckyih, • > • Jk-i). (12) 


and based on this property of cumulants, for the linear and 
time-invariant systems it can be written that 

CO oa 

yd) = ^ x(i - j)h(j) + n(i), ^ \h(i)\ < oo (13) 

j—~CO /=—CO 

and k-th-order cumulants becomes Gil 


Cky(i\,i2, h, ■■ ■, ik-l) - 
oo 

Ckx(0 ,..., 0) ^ h(i)h(i + ii)... h(i + ik-i). (14) 

/= —(Xt 

Assuming that the x(i) is i.i.d., in case of non-Gaussian sig¬ 
nals, zero-lag fourth order cumulants will converge to the finite 
moments of C 4 x{Q, 0,0) 4 0 and in case of Gaussian distributed 
signals C4 j:( 0, 0,0) = 0. Thus C 4 y{Q, 0,0) will also converge to 
zero theoretically. However, this is not the case for wireless 
signals because, the signal model given in Section HD implies 
that wireless channel is time-variant across multiple symbols. 
Moreover, phase and frequency offset can not be assumed 
constant over the period of N symbols which can extend to 
a couple of hundreds as indicated in Section IIII-AI Therefore 
it is not possible to shift the channel coefficients, phase and 
frequency offset components out of the cumulant equations. 
The i.i.d. property will not perfectly hold because of these 
additional components introduced to each sample of the signal 
when compared to eqn.([T]i. Hence the high order cumulants of 
wireless signals will not only be determined by the transmitted 
signal, but also will depend on the communications medium 
characteristics. Furthermore, the probabilistic convergence of 
estimated cumulants Cky to the Cky can be possible when y(i) 
samples are independent and well separated in time. This is 
called mixing conditions. However, the autocorrelation C 2 y(ii) 
is estimated using sample averaging because instead of the 
series y(t), the samples of the original series y(i) are employed 
in the process. Moreover, there is no synchronization process 
involved in the sampling stage. Therefore an irreducible error 
is introduced to the estimation of the cumulants. Under these 
conditions q, will not converge perfectly to zero and can take 
different numerical values based on the channel conditions and 
communications medium. On the other hand, when the fourth 
order cumulants are written in the form of 


C4y(ii, h, h) = E{y(0y(i + ii)y(i + i 2 )yd + h)} 

- E{g(/)g(i ii)g(i + i2)g{i + * 3 )), (15) 

where {g(/)) is a Gaussian random process, it can be inferred 
from the eqn. (m that fourth order cumulants provide a 
measure of distance of the distribution of the signal samples 
form Gaussianity. Considering the fact that the transmitted SC 














5 






2 N-A 


2 ;(A 9 ,+A 9 M)g; 2 ;rA// 


N ■ 


' YfiiAt 





^ L L -|2 

_ ^g2XAe,+ASi+,)g727rA/,(,+i) _ rjhJiAt) ^x((/ + /l)Af - T,„)h,„{{i + A)At) 

i=l m=l m=l -* 


(16) 


signals will also be affected from same propagation character¬ 
istics and consequently their estimates of HOS will change, 
it is still plausible to assume that estimates of the fourth 
order cumulants of OFDM signals will converge to closer 
numerical values to zero when compared to SC carrier signals. 
However, under wireless channel conditions and propagation 
characteristics, the statistical identification methods which 
are based on fixed thresholds or empirical numerical values 
such as 0, 0, HI will fail to identify OFDM signals 
consistently. To this end, a new approach is required to the 
identification of OFDM signals for wireless communications 
systems. 

B. The Gaussianity Test and Decision Mechanism 

It is known from that asymptotic Gaussianity of the 
summable cumulants leads to 

ViV(c4,-C4v) D''(0,s,) (17) 

N-^oo 

where D'' denotes real Gaussian distribution, C 4 y is the 4th- 
order cumulant vector which holds the theoretical cumulant 
values C 4 y A limA^_<co and Ec is the asymptotic covari¬ 

ance matrix of C 4 y which has the close form of 

Sc A lim A^E{(C4y - C4y)(C4y “ C4„)'). (18) 

N->oo . . . y 

where ' denotes the transpose operation. Therefore, the time- 
domain Gaussianity test can be formulated as a binary hypoth¬ 
esis problem of 

-Ko : C4y - D''(0,iV-'Sc) 

VS. 

-Ki : C4y - D''(C4y,JV-%), -- C4y + 0. (19) 

To be able to solve the problem, in this study, the dQ 4 test 
which was proposed in 11211 is adopted to the peculiar case of 
OFDM signal identification under multipath wireless channels. 
Therefore the binary hypothesis problem which is given in eqn. 
(fl^ can be addressed by a chi-square test called as dc ,4 which 
is defined as 

dcA - NC4ytT^^C4y, ( 20 ) 

and from ED it can be seen that 

dcA^m,^t-^C4y X\- ( 21 ) 

The pseudo-inverse of £ ' should replace inverse in the 
equations below if £ rank-deficient. Therefore, for an a level 
significance the test in eqn. (flOl l becomes a chi-squared test 
which can be defined by 

dcA % tc -XN^ia)- ( 22 ) 


tc is disribution table value with degree of freedom of a. 
The probability of the false alarm is given by 

PF^a<P[dcA>x\YHo\ (23) 

Under the distribution of dG,A is estimated from eqn. ([T^ 
as 

dcA ~ D''{Nc\yt^^^C4yANc^^tr^^C4y] (24) 

and the probability of detection is given by 

PD^a<P{dGA>tGYHil (25) 

The fourth order cumulant and covariance matrix estima¬ 

tion processes derived herein is utilized by the introduced 
asymptotic chi-squared CFAR test. Thus the effects of fading 
channels which can lead to stronger non-Gaussian compo¬ 
nents are minimized by the involvement of covariances of 
the cumulants. Selection of the threshold is conducted by 
employing the test tables based on the degree of freedom 
of the distribution. On the other hand, fourth order cumulant 
covariance matrix estimate computation is a long and complex 
procedure which requires involvement of products of the cross 
moment terms. Moreover, computation of the Ec can become 
impractical when C 4 y vector becomes very large. However, 
peculiar to the OFDM signals, under the conditions defined in 
Section ITII-AI covariance matrix estimate calculation becomes 
a simpler process. The general form of covariance estimates 
of the cumulants for OFDM signals is derived based on the 
computations in 1271 (eqn. 2.3.8) and 1211 and reduced to three 
equations which are consist of several moments of the received 
signal. The computation process can be found in the Appendix. 
Furthermore, the test is originally applied to real processes and 
in this paper the same approach is maintained because the 
distributional characteristics of the complex signals are also 
retained by the imaginary and real parts. Besides, operating in 
real or imaginary domain reduces the computation time for the 
test output when compared to complex domain. Therefore, the 
analysis for wireless OFDM signals can be given over either 
of these domains and the analysis can be extended to com¬ 
plex domain l22l . The proposed OFDM signal identification 
method can be executed as follows: 

. calculate the estimates of the fourth order cumulants of 
the received signal based on given degree of freedom, i.e., 
M, using equation (fTTT l. 

. calculate the of the fourth order cumulants based on 
the equations given in Appendix, 

. calculate the chi-squared test result from the equation 
Eii, using the estimates computed in the previous two 
steps, and 
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Parameter 

Value 

Carrier frequency 

2140 MHz 

Tx/Rx distance 

18.24 feet 

Modulation type 

SC-PSK, SC-FSK, SC-QAM, OFDM-QAM 

Modulation order 

2,4,8,16,32,64,128,256 

Number of symbols 

64,128,256,512,1024,2048 

SNR level 

-4 to 15 dB 

Symbol rate 

250 kHz 

Filter type 

Root raised cosine, a = 0,3 

Frequency Offset (A/) 

10 Hz (max.) 


TABLE I 

Measurement setup: the values of the used parameters 


. compare the test result with the tabular value of the test 
with the given degree of freedom and make the decision 
based on equation (l22l l. 

The calculation of the estimates of the cumulants at the 
first step of the proposed method is achieved by a set of mul¬ 
tiplication, shifting, and summation operations for a limited 
number of lags. Moreover in the Appendix, the estimation 
of the covariance matrix in the second step of the proposed 
method is reduced down to calculation of a set of high- 
order moments which does not include any shifting operations. 
These mathematical operations can be implemented over mi¬ 
croprocessors or micro-controllers of the measurement devices 
by the utilization of smart mathematical algorithms which 
avoid calculation duplications. Therefore, the implementation 
of the proposed method can be greatly simplified. Furthermore, 
the construction of the decision process at the last two steps 
consist of very straightforward multiplication and comparison 
operations which do not introduce an additional complexity 
to implementation of the proposed method over the wireless 
communications systems. 

IV. Measurement Setup 

a measurement setup is developed to investigate the identifi¬ 
cation performance of the proposed method. The measurement 
system consists of an Agilent ESG E4438C signal generator 
(Tx), an Agilent E4440A PSA series spectrum analyzer {Rx) 
and accompanying vector signal analyzer (VSA) software to 
record the signals, a tuple of omni-directional antennas to 
cover the wireless communications bands of interest, a laptop 
computer, cables, and connectors. All the measurement pro¬ 
cess is controlled by the laptop computer. The signal generator 
and the spectrum analyzer are connected to the network router 
and the laptop computer commands both devices through 
the transmission control protocol and internet protocol. The 
signal transmission and the reception starts and ends with 
the commands that are sent by the laptop computer. The 
setup is depicted in Eig. |l(a)[ Moreover, the signal generator 
supports the generation and transmission of waveforms that 
are designed separate from the generator. Therefore, a signal 
generation software is developed in MATLAB to be able to 
obtain the full control of the data generation and transmission 
procedures. 

The developed software generates the SC-PSK, SC-ESK, 
SC-QAM, and OEDM modulated signals with various mod¬ 
ulation orders based on the signal model in Section |II] The 


OEDM signals employ QAM as the underlaying modulation 
technique over 64 active catTiers with 250 kHz spacing. Tq 
is also selected as IjATr,. The signal generation software 
generates the designed waveforms at the laptop and uses ven¬ 
dor provided library functions to access the signal generator. 
As the software accesses the generator, it uploads generated 
waveforms to the device and commands the signal generator 
to transmit the predefined number of symbols of the waveform 
at the defined frequency band, with the predefined signal 
power. Moreover, signal parameters such as modulation type 
and order can also be adjusted by the software. Therefore 
the parametric analysis which is based on the measurement 
parameters such as SNR, modulation type, data length, and 
modulation order can be done methodologically. On the other 
hand, the VS A software provided by the vendor is employed 
for signal acquisition and recording purposes. SNR of the 
received signal is also determined by this software. Therefore, 
all the signal transmission loop is controlled and completed 
on the laptop computer while the transmitter and the receiver 
is allocated to different locations in the laboratory. The list of 
the all parameters used in the measurement process is given 
in TABLE U The measurements are repeated for 200 times 
for each changing parameter by the aid of signal generation 
software. 

The measurement environment also has an important effect 
over the performance of the proposed method. The TxjRx 
devices are fixed to certain allocations in the measurement 
laboratory which has scattering surfaces such as metal panels, 
other measurement and calibration equipment, tables, chairs 
etc., as shown in Eig. |l(b)| Therefore, the measurement envi¬ 
ronment is a highly reflective indoor environment. It should be 
noted that the performance of any signal identification method 
will be affected under these conditions, because the wireless 
channel is defined by the heavily reflective nature of the indoor 
environment. These circumstances can be assumed to be the 
bottom line performance conditions and it can also be assumed 
that the identification performance of the proposed method will 
increase as the system is deployed outdoors with relatively less 
reflections. 

V. Measurement Results 

a gear detector is implemented based on the proposed 
method and the equations (l20l) - (l25]) in Section Ull-Bl for Pf of 
0.001, 0.01, 0.1 and 200 measurements. On the other hand, a 
frequency shift of A/ = 10 Hz is added to the signal generation 
process to imitate the Doppler shift and other impairments that 
cause the shift in the frequency. In Eig. |2] [3] and |4] the Pa 
results for OEDM vs. SC signals are given for 1024 symbols, 
modulation order of 32, M - l.STs, and various SNR levels. 
The results indicate that best separation is achieved against 
SC-ESK, SC-PSK and SC-QAM consecutively. Especially for 
Pf of 0.1, and 0.01, a certain confidence level is achieved at 
2 dB of SNR for SC-ESK signals, at 6 dB for SC-PSK, and 
at 9 dB for SC-QAM. 

The detection performance of the proposed method falls 
sharply in a range of couple of dBs at the end of positive 
SNR region for all signals and identification becomes impos¬ 
sible when the SNR is negative. The measurement results 
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Vector Signal 

Signal Generator Analyzer 



(a) Measurement setup: transmitter, receiver, laptop com¬ 
puter, router and their connections. 


Fig. 1. The measurement setup and environment. 



(b) A general overview of the measurement environment. 



Fig. 2. Probability of detection for SC-FSK v.s. OFDM: N = 1024, 
modulation order = 32, M = l.5Ts. 



SNR [dB] 


Fig. 4. Probability of detection for SC-QAM v.s. OFDM: N = 1024, 
modulation order = 32, M = l.5Ts. 



SNR [dB] 


Fig. 3. Probability of detection for SC-PSK v.s. OFDM: N = 1024, 
modulation order = 32, M = l.5Ts. 


indicate that the noise dominates the test statistics in this 
region, therefore identification performance deteriorates. SC- 
QAM follows other modulation types in terms of performance 
in all measurements because the underlying modulation that 


is employed by the OFDM signal is also QAM. For the 
modulation order of 32, 1024 symbols, and Pf =0.1, detection 
performance above 0.9 is achieved at 0 dB of SNR for OFDM 
vs. SC-FSK, at 2 dB for OFDM vs. SC-PSK, and at 3 dB for 
OFDM vs. SC-QAM. 

The measurement results are compared with the simulation 
results which are based on ll23l for each signal type and 
Pf — 0.01 in Fig.|5]|6] and|7] Except for a set of negative SNR 
values, in general the comparison plots indicate that simulation 
results are slightly better and in some cases -especially in 
the initial positive SNR region, an increase of 0.1 at the Pq 
values can be observed for all signals. This is the indication 
of the significant effect of measurement environment depicted 
in Section |IV] over the wireless channel and consequently 
on the performance of the proposed method in contrary to 
the four taps simulated channel implemented in Il23l . On the 
other hand, transmitter and the receiver are directly connected 
with a cable and the measurement process is synchronized by 
the available software libraries in both of the devices. These 
conditions imply only the effect of AWGN is introduced to 
the received signal during the transmission, thus the plots 
that are given in the same figures with the label of AWGN 






























































































SNR 


Fig. 5. Pd performance comparison for SC-FSK v.s. OFDM: AWGN channel, 
simulations, measurements. {Pf = 0.01), N = 1024, modulation order = 32, 
M = \.5Ts 



SNR 


Fig. 6. Pd performance comparison for SC-PSK v.s. OFDM: AWGN channel, 
simulations, measurements. {Pf = 0.01), N = 1024, modulation order = 32, 
M = l.5Ts 


can be considered as the upper bound performance limits for 
the proposed signal identification method. These measurement 
indicate Pq above 0.9 for 0 dB of SNR in case of FSK vs. 
OFDM, 1 dB in case of PSK vs. OFDM, and 5 dB in case of 
QAM vs. OFDM. 

The Pd performance of the proposed method is also in¬ 
vestigated by employing QPSK modulated SC wideband code 
division multiple access (WCDMA) and OFDM based LTE 
signals. The transmission bandwidth for WCDMA is 5 MHz 
and LTE system is also conhgured to 5 MHz transmission 
bandwidth with 15 KHz spaced carriers and fast Eourier trans¬ 
form size of 512. The signals are generated from the wireless 
standards based libraries available in the Agilent ESG E4438C 
signal generator. The rest of the measurement parameters are 
kept same as given in TABLE|I1 The measurement results show 
similarity with that of generic SC-PSK measurements however, 
Pd leads to certain identihcation at higher SNR values when 
compared to the generic case. 

The effect of the change of modulation order on the test 



Fig. 7. Pd performance comparison for SC-QAM v.s. OFDM: AWGN 
channel, simulations, measurements. (Pf = 0.01), N = 1024, modulation 
order = 32, M = l.5Ts 



SNR [dB] 


Fig. 8. Probability of detection for WCDMA v.s. LTE: M = l.STg. 


statistics is investigated in Pig. |9] when the SNR is 7 dB, 
Pf - 0.1, and 1024 symbols are recorded. While SC-PSK 
separation barely depends on the modulation order, as the 
order increase above 8, the proposed method performs better 
for SC-PSK steadily, and Pd reaches to the levels above 0.95. 
In case of SC-QAM signals, the best performance is achieved 
as the modulation order reaches to 32 with Pd > 0.96. 
The measurement results indicate no significant effect of 
modulation order over identification performance. 

Data length or number of symbols used is also another 
important parameter which affects the performance of the 
statistical tests. Pig. [TO] displays the detection performance of 
the proposed method with the changing number of symbols 
for 7 dB SNR, P/r = 0.1, and modulation order of 32 for each 
signal type. Below 128 symbols the proposed method performs 
poorly due to the problems discussed in Section |III] While 
the test performs well for PSK signals over 256 symbols with 
Pd > 0.95, for all SC signals the test performance is optimized 
as the number of symbols reaches 512. The degradation in the 
test performance after 1024 symbols is due to central limit 












































































































































































































9 



Fig. 9. Effect of modulation order on the identification: {Pf =0.1), SNR = 
7 dB, W= 1024, M= \.5Ts. 



Fig. 10. Effect of number of symbols on the identification: (Pf = 0.1), SNR 
= 7 dB, modulation order = 32. M = 1.57^. 

theorem. 

The lag distance which also determine the degree of free¬ 
dom of the test has an important effect on the identification 
performance of the proposed method. More information about 
the distribution of the received signal is introduced to the 
decision process as the lag distance increases. However, the 
computation time required to calculate the estimates of the 
cumulants will also increase as the lag distance increases. 
The detection improvements for increased lag distance in 
terms of symbol duration is given in Fig. [TT]for N = 1024, 
SNR = 7 dB, modulation order of 32 and Pp - 0.01. The 
oversampling rate is double the Nyquist rate, therefore, for 
instance, M - Q.5Ts corresponds to 2 samples. A significant 
gain of 0.4 in the detection performance can be noticed in 
case of QAM signals as the lag distance increase from 0.5 
to 2.5, however computation time requirements also increase 
significantly. For other signal types the Pd performances are 
steady Le., above 0.9. 

Finally the identification performance of the proposed 
method is compared with three of the non-coherent statistical 
methods and a coherent method that are available in the 



0.5 1 1.5 2 2.5 

Lag Distance (M) [Symbol duration (Tg)] 

Fig. 11. Effect of lag distance (degree of freedom) on the identification: 
{Pf = 0.01), N = 1024, SNR = 7 dB, modulation order = 32. 


literature. One of these methods, which is introduced in la 
depends on the fourth and the sixth orders of moments of 
the OFDM signals and introduces an MMSE algorithm based 
on these moment functions (MMSE-M). Secondly in la, a 
method based on the combination of fourth order moments 
and the amplitude estimation is introduced (AE-M) and third, 
sums of fourth order cumulants are used for OEDM v.s. SC 
classification in ||4| (SUM-C). The coherent method employs 
a two stage identification procedure to distinguish SLCD, 
BT-SLCD, and OEDM signals from each other based on 
the unique second order cyclostationary features of these 
signals. Eirst the SLCD signals are distinguished from the 
BT-SLCD/OEDM pair, then these two signals are identified 
in the second stage HSl (BT-SLCD), however the first stage 
of the signal classification algorithm is investigated for the 
comparison purposes assuming that the there are two possible 
signal types. The performance comparison with these methods 
can be found in Fig.[T2]for the SC-PSK case. Modulation order 
is selected as 32, Pp - 0.01, and N - 1204. 

The proposed method performs better than the non-coherent 
methods because of the underlying test structure which does 
not depend on fixed thresholds and utilization of the estimates 
of the covariances of the cumulants in the identification pro¬ 
cess. At 1 dB SNR region, the proposed method outperforms 
other methods with the Pd margin of 0.2. Eor the 5 dB of 
SNR, the margin is 0.15 with SUM-C, 0.08 with AE-M, and 
0.05 with MMSE-M. As the SNR reaches to 10 dB, the 
tests converge to closer Pd rates. On the other hand, while 
the coherent BT-SLCD algorithm outperforms the proposed 
method on the low-SNR region, the identification performance 
of the proposed method converges to the BT-SLCD around 
the 8 dB SNR region. These results are expected when it 
is considered that in general coherent techniques have better 
performance than the statistical test, trading a certain level of 
increase in the algorithm complexity and availability of certain 
information to the receiver beforehand with improvements in 
terms identification performance. The performance compari¬ 
son for SC-ESK and SC-QAM signals are not included, but 
similar results are observed. 
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Fig. 12. Performance comparison with other algorithms: SC-PSK v.s. OFDM, 
{Pf = 0.01), N = 1024, modulation order = 32, M = \.5Ts. 


cumulants can be written in terms of moments as ED 


C4y(ii,i2, h) = myiiuh, h) - m2y(ii)m2y(i2 - h) 

- m2y(i2)m2y(h - h) - my(h)m2y(ii - h) 

(A.l) 


where 

myik) - + 'iMf + ' 2 )■ ■ -yit + ik-\), (a.2) 

and 4 = (ii,i 2 ,. ■. ,ik-i)- The computation of each entry 
cov{c 4 y(ii,i 2 ,h),C 4 y(ji,j 2 ,j 3 )] of the 2 c becomes a compu¬ 
tationally consuming process under . However the compu¬ 
tation of the diagonal entries defined by is sufficient for 
OFDM signals. If d = ii - ( 2 , P - ji - 72 . h - 73 = 0 and 
0 < P < M, the open for of covariance of Qy can be given in 
three parts and the first part is consists of covariances of two 
fourth order moments as 


VI. Conclusion 

Signal identification is a vital component of modern wireless 
communications systems. In this paper first, fourth order 
cumulant estimates of wireless OFDM signals are given. Sec¬ 
ondly, it is shown that wireless impairments have significant 
impacts on the distribution characteristics of wireless signals. 
Then, an OFDM signal identification method which employs 
estimates of fourth order cumulants and their covariances for 
limited number of lags is proposed. The cumulant and the co- 
variance matrix calculations are greatly simplified particularly 
for OFDM signals. The method utilizes a modified test at 
the decision process and symbol duration information is used 
to determine the number of cumulant legs or the degree of 
freedom of the underlying Gaussianity test. In practice, degree 
of freedom which is M > 10 can satisfy certain confidence 
levels, thus the requirement for the symbol duration can also 
be waived. A measurement setup is developed to validate 
the performance of the proposed method realistically. The 
identification performance of the method is provided utilizing 
parameters such as modulation order and type, data length, 
and degree of freedom of the test. A CFAR detector is 
used in the performance analysis. The method performs better 
than other statistical identification methods under the wireless 
multipath fading channel conditions and it is robust against 
the changing modulation order and data length. However the 
identification performance is affected under low-SNR, because 
of the time-domain Gaussianity of the OFDM signals. Pd of 
more than 0.9 is achieved for all signal types when Pf =0.1. 


Appendix 

The asymptotic covariances of sample cumulants can be 
estimated based on the straightforward asymptotic covariance 
expression for the cumulant orders up to A: = 3 Il2^ . However, 
the covariance estimation becomes a complex procedure for 
k > 3 due to the combinations of products of the second 
and fourth order moments. The estimate of the 4-th order 


Cl - cov{m4y(A, A, 0 ),m 4 y(p,p, 0 )}, (A.3a) 

and the second part will be consist of various combinations of 
covariances of fourth order and second order moments. The 


open expression of these terms are given by 

C 2 - cov{m4y(A, A,0), m 2 y(P)m 2 y(P)} (A.3bl) 

+ cov{m4y{A,A,0),m2y(P)m2y{-p)} (A.3b2) 

+ cov{m4y{A,A,Q),m2y{G)m2y(Qi)} (A.3b3) 

+ cov\m4y(p,p,0),m2y{A)m2y{A)} (A.3b4) 

+ cov{m4y{p,p,0),m2y{A)m2y{-A)} (A.3b5) 

+ cov{m4y{p,p,Q),m2y{Q)m2y{Q)}, (A.3b6) 


and finally the third term will be constituted by the covari¬ 
ances of all possible second order moments of the estimated 


cumulants as 

C 3 - cov{m 2 y(A)m 2 y(A),m 2 y(P)m 2 y(P)} (A.3cl) 

+ cov{m 2 y(A)m 2 y(A),m 2 y(P)m 2 y(-P)} (A.3c2) 

+ cov{m2y(A)m2y(A),m2y(0)m2y(0)} (A.3c3) 

+ cov{m2y(A)th2y(-A), m2y(P)m2y(P)] (A.3c4) 

+ cov{m2y(A)m2y(-A), m2y(P)m2y{-p)] (A.3c5) 

+ cov{m 2 y(A)m 2 y(-A), m2y(0)m2y(0)} (A.3c6) 

+ cov{m 2 y(Q)m 2 yiQ), m 2 y(P)m 2 y(P)} (A.3c7) 

+ cov{m2y(0)m2y(0), m 2 y{p)m 2 y{-P)} (A.3c8) 

+ cov{m 2 y{Q)m 2 y{Q), m 2 y(Q)m 2 y{Q)]. (A.3c9) 

Therefore the covariance of C 4 y is given by 

cov{c4y(A, A, 0), C4y(p,p, 0)) = Cl -I- C2 -I- C3. (A.3d) 


The derivation of these three groups can be based on 
the computable sample approximations of the asymptotic 
covariances established in ll3TII . Il^ for the mixing discrete¬ 
time stationary processes such as defined in eqn. (fT4l i. The 
covariance of the moments is defined as 

c» 

lim Ncovimydic), mkyi ji)] = / , < 2 ov{y( 0 )y(ii).. .y(4-i), 

A'—>00 ' / i 

—00 

y(T)y(T + ji )... y(T + ji - i )}. (A.4) 
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Following cov{ab} = E{ab} - E{a}E{b} based on eqn. ( IA.4l i 
and replacing the moments with their estimates as defined in 
eqn. ( IA.2I ). (A.3a) can be written as 

Cl = cov{m 4 y(A, A,0),m4y(J3,/3,0)} 

~ ^ [m&y(A, A,0 ,t,t+/3,t+ jA, t ) 

t=-Kn 

- OT4y(/l, d, 0)m4v(y6,jS, 0)] (A.5) 

Limits of the sum in the eqn. (IA.5I) should satisfy 
as N ^ oa and K^fN —> 0 as A —> oo. In practice, for OFDM 
systems Kf^ can be selected as Kf^ > t where C 2 .i„,(t) ~ 0 ||2TI| . 
The second group of sample estimates include the covariances 
of the multiplications of the moment terms and theorem 2.3.2 
in Il27ll expresses such terms in the form of 

cov{a, be] - cov{b}cov{a, c) + cov{c}cov{a, b] (A.6) 

and eqn. (IA.3ab can be written as 

cov{m4y(d, A, 0), m2y(Ji)m2y(Ji)) 

- E{m2y{p)}cov\m4y{A, A, 0), m2y(j3)} 

+E{m2y(j8)}cov{m4y(A, A, 0), myiJA)} 

- 2^E{m2y(JA)}cov{m4y(A, A, 0), myiJA)}]. 

(A.7) 

Depending on the estimator given in eqn. (IA.21 i and covariance 
definition in eqn. (IA.41 i. similar to eqn. (IA.51 l. the covariances 
in eqn. (IA.7b can be expressed in terms of moments as 

cov{m4y(A, A, 0), m2y(JA)m2y(/A)} 

= 2[m2y(J3) ^ [mey(A, A, 0,t,t+ (A)- myiA, A, 0)m2y(JA)]]. 

r=-Kn 

(A.8) 

When the same derivations are applied to the rest of the 
equations and after some simplifications (IA.3al - IA.3b6b can 
be written as 

C2 = 2{p.{2m2y{JA) + m2yi-/A)) + <^m2y(JA) + 2'I'OT2y(0)), (A.9) 
where 

— ^ [mey(A,A,0,T,T+P)-m4y(A,A,0)m2y(JA)], (A.IO) 
t=-Kn 

O = — ^ [m6y(d, d, 0, T, T - yS) - OT 4 v(d, d, 0)m2y{-/A)], 
t=-Kn 

(A.ll) 

and 

'T = — ^ [m6,,(d,d,0, T,T)-m4y(d,d,0)TO2y(0)]. (A. 12) 

t=-Kn 

The leads to symmetric Sc therefore for each Sc entry not 
all six terms of second groups should be processed, estimation 
of the first half would be sufficient and the following terms will 
be replicated as indicated in eqn. (IA.9b . Similar to the second 
group, covariances of multiplications of the estimate moments 
which consist the third group are in the form of cov{ab,cd}. 


According to eqn.(2.3.8) in IIZTll these terms can be expressed 
as 

cov{ab,cd} - cov{a,c}cov{b,d} + cov{a,d}cov{b,c}, (A. 13) 

then the eqn. (IA.3b6b is given by 

cov{m2y(A)m2y(A),m2y(JA)m2y(JA)} 

= cov{OT2y(4), m2y(fi)}cov{m2y(A), m2y{jA)] 

+ cov{m 2 y(A), m 2 y(JA)}cov{m 2 y(A), myifA)}. (A. 14) 

After applying the conversion defined in the eqn. (IA.13b and 
rearranging all nine terms, (IA.3b6l - IA.3c9b can be written as 

C 3 =2r^ + 2FA + 20^ + 2Fcov{m2y(d), rh 2 yi-/A)} 

+ A(r + A) + 2©cov{m2y(-d), m2y(0)) + 2r^ 

+ 2Tcov{m2y(0),m2y{-fA)} + 2[cov\m2y{0),m2y{Q)}f, 

(A. 15) 

where 

F=— [m4y(A, T, T + /?) - m2y(A)m2y(fA)] , (A.16) 

t=-Kn 

A = — ^ [m4y(A, T,T- /A)- m2y(A)m2y(-/A)], (A. 17) 

t~-Kn 

1 

©=— y [myjA, T, t) - m2y(d)m2y(0)], (A.18) 

^ r=-Kf, 

A = — ^ [m4y(-A, T, T -jS) - m2y(-A)m2y{-/A)], (A. 19) 

^ Kn 

T = — ^ [myiA, T, T + y6) - OT2y(4)m2y(y6)], (A.20) 

^ t=-Kh 

and 

COV{OT2y(-d), myiJA)} 

2 ^A 

= — 'y[m4yi-A,T,T+IA)-m2yi-A)m2yi[A)], (A.21) 

^ T=-Ky 

cov{m2yi-A),m2yiO)} ^ ^ [m4y(-d, r, T)-m2y(-d)m2y(0)], 

T=-^A' 

(A.22) 

cov{m2yiO),m2yi-/A)} ^ ^ [m4y(0, r, T-/3)-m2y(0)m2y(-j8)], 

^ t=-Kh 

(A.23) 

j ^A 

COv{m2y(0), m2y(0)) = — y [OT4y(0,T,T) - OT2y(0)OT2y(0)]. 
t=-Kn 

(A.24) 

Each entry of Ec now can be computed as the sum of equations 
(IA.5b . ( IA.9b . and (IA.15b under 
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